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Abstract 

We present a systematic method to determine BCJ numerators for one-loop amplitudes that 
O explores the global constraints on the loop momentum dependence. We apply this method to 

amplitudes in = 4 gauge theory, working out detailed examples up to seven points in both the 
MHV and the NMHV sectors. The structure of Jacobi identities between BCJ numerators is seen to 
^ be closely connected to that of algebraic integrand reductions. We discuss the consequences for one- 

loop J\f = 8 supergravity amplitudes obtained through the double copy prescription. Moreover, 
^ in the MHV sector, we show how to obtain simple BCJ box numerators using a conjectured 

relationship to amplitudes in self-dual gauge theory. We also introduce simpler trace-type formulas 
for integrand reductions. 
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I. INTRODUCTION 

In 2008, Bern, Carrasco, and Johansson uncovered a curious property of tree- level am- 
plitudes in A/" = 4 supersymmetric Yang-Mills (SYM) theory [1]. They noticed that it is 
possible to rearrange the amplitudes into a form where kinematic factors and colour factors 
appear symmetrically, and where the kinematic factors satisfy the same algebraic identities 
as the colour factors. This property has come to be known as the BCJ duality between 
colour and kinematics, and it suggests the presence of a hidden symmetry in a wide range 
of Yang-Mills theories. Moreover, one can construct certain gravity amplitudes by replac- 
ing colour factors with another copy of BCJ dual kinematic factors, exposing a mysterious 
connection between the structure of amplitudes in Yang-Mills and gravity theories. 
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Since their paper came out, many interesting features of the tree-level duality have been 
found [2-13]. The so-called BCJ amplitude relations, new linear relations among tree-level 
partial amplitudes that follow from the duality, add to the previously known Kleiss-Kuijf 
relations [14] and have been proven using a variety of methods [15-19]. The double copy 
relation to gravity has also been proven at tree level [5]; it is equivalent to the celebrated 
KLT relation [20], proven in field theory for the first time in [21]; see also [22, 23]. 

However, what is most striking about the BCJ duality in Yang-Mills theory and its double 
copy relation to gravity is that they seem to hold also at loop level [24]. They have proven 
to be a very useful tool for studying multi-loop amplitudes in jV = 8 supergravity [24-27], 
A/" = 4 supergravity [28-31], and also theories with less supersymmetry [32]. Additional 
evidence for this loop-level structure appears in soft [33] and high-energy [34] limits, and 
also in the BCFW shifts of gauge theory integrands [35]. Relations between gauge theory 
and gravity at one loop were also investigated in [36], and a one-loop version of the BCJ 
amplitude relations has been proven in [37]. 

Because of the connection between kinematics and colour inherent in the BCJ duality, 
one might expect an underlying kinematic symmetry to be at play. Indeed, a partial un- 
derstanding of the kinematic symmetry has appeared in refs. [9, 12], and this was later 
used to construct a BCJ representation for finite one-loop amplitudes in pure Yang-Mills 
theory [38]. Beyond Yang- Mills scattering amplitudes, the duality has also appeared in 
multi-loop form factors [39] and higher-dimension operators [40]. There is even evidence of 
the duality holding in ABJM theory, where the colour is associated to a 3-algebra instead 
of a (Lie) 2-algebra [41, 42]. 

One challenge to further progress in the gauge theory amplitude context is that it is very 
difficult to obtain BCJ dual forms of amplitudes at high loop orders and high multiplicities. 
There is no known construction in general; instead, the usual strategy is to start with an 
ansatz for a set of 'master' numerators, from which all other numerators can be obtained 
through Jacobi identities. One then fixes unknown coefficients by matching the ansatz to 
unitarity cuts of the amplitude [43-46]. The ansatz grows combinatorially as the loop order 
and multiplicity are increased, which puts a practical limit on the problems this can be 
applied to. Moreover, it is always possible that the ansatz is not general enough to contain 
the actual BCJ form of the amplitude. 

For these reasons, we develop here a new strategy for obtaining BCJ numerators con- 
structively rather than by ansatz. Our strategy, which we apply to A/" = 4 SYM amplitudes, 
utilizes global properties of the system of one-loop numerator Jacobi identities to invert the 
usual top-down approach; instead of starting with master numerators and constructing all 
other numerators via Jacobi identities, we start with the much simpler box numerators and 
construct all other numerators in a bottom-up manner, a procedure we term "integrand 
oxidation." This still assumes knowledge of the BCJ box numerators, and we present a 
method for constructing these in the case of MHV amplitudes. Moreover, in the absence 
of such BCJ box numerators, for example for non-MHV amplitudes, we present a method 
starting from unitarity cuts of the amplitude. 

This paper is organised as follows: In Section II, we review the BCJ duality between colour 
and kinematics. In Section III, we tackle the case of the five-point one-loop MHV amplitude 
in A/" = 4 SYM theory. Although this has been thoroughly discussed elsewhere [26, 47, 48], 
we take the opportunity to introduce our general strategy with the simplest possible ex- 
ample. We see explicitly that knowledge of BCJ boxes in this case is enough to fix the 
entire amplitude. In Sections IV and V, we discuss two ingredients that will be necessary to 
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progress to higher-point amphtudes. First, we give an exphcit construction of BCJ box nu- 
merators by exploiting a dimension-shifting relationship between one-loop MHV amplitudes 
in A/^ = 4 SYM theory and one-loop all-plus amplitudes in pure Yang-Mills [45]. Then we 
discuss scalar integral reduction techniques, presenting simplified expressions. In Section VI, 
we use the six-point box numerators obtained in Section IV to reconstruct the pentagon and 
hexagon BCJ numerators, including all loop momentum dependence. We build on this result 
in Section VII by constructing the BCJ form of the seven-point MHV amplitude in a similar 
manner before progressing to non-MHV six- and seven-point amplitudes in Section VIII. We 
wrap up with a discussion of how cancellations in gravity amplitudes arise from the BCJ 
form of the amplitude in Section IX, and present our conclusions in Section X. 



II. REVIEW 

The statement of the BCJ colour-kinematics duality is that gauge theory amplitudes can 
always be written as a sum over cubic Feynman-like diagrams, whose kinematic dependence 
mirrors the algebraic properties of the colour factors. Consider tree-level amplitudes. Associ- 
ated with each cubic diagram Fj are three objects: a colour factor q, a kinematic numerator 
rij, and a set of propagator denominators P^i where k runs over the set of propagators in Fj. 
The amplitude is 

A^ = g-'^ E (2.1) 

diagrams Ti ^ ' 

Comparing with the usual Feynman diagram expansion in gauge theory, one can see that 
to arrive at Eq. (2.1), one must absorb diagrams involving four-point vertices into the nu- 
merators rii. This can always be achieved by multiplying and dividing contributions by a 
propagator factor. This is not a unique procedure, so there are many sets of numerators 
such that Eq. (2.1) is a valid expression for the amplitude. The key discovery of BCJ was 
that certain representations originating from such reshufflings are special. They involve 
numerators rii satisfying 

rii ± rij ± Uk = , (2.2) 
whenever the associated colour factors satisfy a Jacobi identity 

Ci±Cj±Ck = 0. (2.3) 

A set of numerators rii which satisfy all the Jacobi relations is said to be colour-kinematics 
dual; for brevity, we will also call them simply dual numerators. 

It is a remarkable fact [1, 5] that replacing the colour factors Cj in Eq. (2.1) with a set 
fii of kinematic numerators leads to something useful as well, namely a representation of 
gravity directly factorized as a double copy of gauge theories: 

M^ = ^{K/2r'^ Y: (2.4) 

diagrams Li -^"^ 

The second set of numerators, fii in this double-copy equation need not correspond to the 
same gauge theory. For example, one can choose rii to be a set of numerators appropriate 
for A/" = 4 super- Yang-Mills, while n, are numerators for non-supersymmetric gauge theory. 
The resulting gravitational theory in this example is A/" = 4 supergravity. The relationships 
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among pairs of gauge theories and a gravity theory estabhshed by the double-copy have 
recently been explored [32, 49]. 

The story we have discussed so far was at tree level. The colour-kinematics duality 
and the double copy are conjectured to extend to all loops [24], but beyond tree level less 
is known. At L loop level, we still express m-point amplitudes as a sum over cubic 
diagrams, 

diagrams Vi j = l ^ ' ^ ^'^ ^ 

where Si is a symmetry factor. Just as at tree level, the condition for colour-kinematics 
duality is 

Ci±Cj±Ck = 0^ ni{i) ± nj{i) ± Ukii) = . (2.6) 

The loop momenta in each numerator are chosen so that the four legs involved in the Jacobi 
identity have the same momentum in the three diagrams. Gravity amplitudes are then 
expressed as 



diagrams 



There is, as yet, no proof that dual numerators always exist at loop level for any theory. 
However, numerators have been discovered in examples. In the case of the four-particle 
amplitude in = 4 gauge theory, dual numerators exist for up to four loops [24, 25] . The 
five point amplitude also exists at up to two loops in the same theory [26]. In the context of 
studying the interplay between the double-copy construction and orbifolds in supergravity, 
the one-loop four point amplitude of both M = 1 and Af = 2 gauge theory were recently 
calculated [32]. The first examples of infinite families of one-loop amplitudes in colour- 
kinematics dual form were discussed in [38]. The relevant families of amplitudes are the 
helicity-equal amplitudes and amplitudes which have m — 1 particles of one helicity, and 1 
particle of the other, in pure gauge theory. It is also known that, at two loops, the four 
particle helicity-equal amplitude can be put in a dual form [24]. 

Our focus in this paper will be on constructing colour-dual numerators for one-loop 
amplitudes in maximally supersymmetric gauge theory and gravity. We now turn to this 
topic in the simplest non-trivial case, namely the five-point amplitude. 



III. AN INVITATION: THE FIVE-POINT AMPLITUDE 

A set of dual numerators at five points was first found by Carrasco and Johansson [26] 
and was more recently discussed by Yuan [47]. We will re-derive this set of numerators in a 
manner that will extend to higher-point amplitudes. 

Suppose we know the BCJ box numerators for this amplitude; for example, we can com- 
pute box cuts in self-dual Yang-Mills theory and extract candidate MHV box numerators, 
as will be described in Section IV. We will label the numerator of the pentagon with legs 
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ordering (a, b, c, d, e) as 725(0, b, c, d, e); then the BCJ box numerators are given by^ 

n,{[l, 2], 3, 4, 5) = n5(l, 2, 3, 4, 5) - n,{2, 1, 3, 4, 5) (3.1) 



^5'(Q) 



[12]2[34][45][53] 



(3.2) 



tr5(1234) 

where tr5(1234) is the usual trace over gamma matrices, tr5(1234) = tr (75^*^1*21*31*4)' 
S^{Q) is the supersymmetric deha function, which gives rise to the gluon MHV factor (ij)^, 
i and j being the negative hehcity gluons. This set of box numerators has the important 
properties that it does not depend on loop momentum, and that it is symmetric for permu- 
tations of its corners, which implies the vanishing of triangle and bubble numerators. 

The Carrasco- Johansson set of numerators is independent of loop momentum. We do 
not expect the same simplification to generalise to higher points (and we will confirm this 
expectation below). Let us therefore open with a set of numerators at five points which may 
depend linearly on the loop momentum i: 



n,{l, 2, 3, 4, 5; i) = ^,5(1, 2, 3, 4, 5) ■ £ + ^0,5(1, 2, 3, 4, 5), 



(3.3) 



where we choose our loop momentum to be between the last and first argument of the 
numerator. That is. 



n5(l,2,3,4,5;£) O 




(3.4) 



The basis of our method is to explore the constraints on the loop momentum dependence 
coming from the Jacobi identities. Consider the following set of equations, obtained by using 
Eq. (3.1) iteratively to move leg 1 around the loop: 



- 715(2, 3, 1,4, 5; £) 

- 715(2, 3, 4,1, 5; £) 
-7i5(2,3,4,5,l;£) 



714(2, [1,3], 4, 5) =715(2, 1,3, 4,5; i 
714(2,3, [1,4],5) =715(2, 3, 1,4,5; i 

714(2, 3, 4, [1, 5]) = 715(2, 3, 4, 1, 5; i) - 715(2, 3, 4, 5, 1; i) . (3.5) 

Summing these equations, we learn that 

n4([l, 2], 3, 4, 5) + 114(2, [1, 3], 4, 5) + 714(2, 3, [1, 4], 5) + 114(2, 3, 4, [1, 5]) 

= 715(1, 2, 3, 4, 5; £) - 715(2, 3, 4, 5, 1; i) . (3.6) 

However, the pentagon numerator must satisfy the cyclic permutation property, 

3 2 



115(2,3,4,5,1;. 





^n5(l,2,3,4,5;£-pi). (3.7) 



^ The subscript on indicates that we are talking about the BCJ numerator of an m-gon, the precise 
identity of which can be read off from the pattern of square brackets surrounding the arguments of nm{- ■ ■)■ 
Although the subscript is redundant notation, we leave it for clarity; however, we stress that the use of n^, 
for example, may indicate different quantities depending on the particular amplitude under discussion. 
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It follows that 



7^5(1, 2, 3, 4, 5; £) - n^{2, 3, 4, 5, 1; = ^1,5(1, 2, 3, 4, 5) ■ pi . (3.8) 

Using this observation in Eq. (3.6), we discover a simple expression for the linear pentagon 
numerator, ni 5: 

n4([l, 2] , 3, 4, 5) + n4(2, [1, 3], 4, 5) + n^{2, 3, [1, 4], 5) + n^{2, 3, 4, [1, 5]) = ^1,5(1, 2, 3, 4, 5) ■ pi . 

(3.9) 

Now, using the Carrasco- Johansson box numerator, Eq. (3.2), we find that the left-hand-side 
of this expression vanishes. This is consistent with the fact that the Carrasco- Johansson 
numerators are independent of loop momentum, so that n5(l, 2, 3, 4, 5; tj = no,5(l, 2, 3, 4, 5). 
At higher points, we will see below that dual numerators do depend on loop momentum. 
Moreover, relations of this type, which fix the tensor numerators in terms of the scalar 
numerators, generalise straightforwardly. Thus, knowledge of the scalar numerators is suffi- 
cient to recover the full i dependence of the numerators. This is a key step in our oxidation 
method. 

Another useful equation can be found in a similar manner, by using the reflection sym- 
metry of the pentagon, 

n5(l,2,3,4,5) = -n5(5,4,3,2,l), (3.10) 

where we no longer consider loop momentum dependence. The sign follows from the an- 
tisymmetric cubic vertices. As before, we begin by moving leg 1 until it sits beside leg 
5: 

n,{[l, 2], 3, 4, 5) = ^5(1, 2, 3, 4, 5) - n,{2, 1, 3, 4, 5) , 
n,{2, [1, 3], 4, 5) = n,{2, 1, 3, 4, 5) - n,{2, 3, 1, 4, 5) , 
n4(2, 3, [1, 4], 5) = ^5(2, 3, 1, 4, 5) - 715(2, 3, 4, 1, 5) . (3.11) 

Our goal now is to move the other legs until we reverse the order of all the legs: 

n4([2, 3], 4, 1, 5) = n,{2, 3, 4, 1, 5) - n5(3, 2, 4, 1, 5) , 
714(3, [2, 4], 1, 5) = ^5(3, 2, 4, 1, 5) - ^5(3, 4, 2, 1, 5) , 
n,{[3, 4], 2, 1, 5) = ^5(3, 4, 2, 1, 5) - n,{A, 3, 2, 1, 5) 

= n5(3,4,2,l,5)-n5(5,4,3,2,l). (3.12) 

The summation telescopes to 

2n5(l, 2, 3, 4, 5) = 715(1, 2, 3, 4, 5) - 715(5, 4, 3, 2, 1) 

= 7i4([l, 2], 3, 4, 5) + 114(2, [1, 3], 4, 5) + 714(2, 3, [1,4], 5) 
+ 7i4([2, 3], 4, 1, 5) + 714(3, [2, 4], 1, 5) + n4([3, 4], 2, 1, 5) 

tr5(1234) ^ ' 

Thus, knowledge of the BCJ box numerators is sufficient to reconstruct the BCJ pentagon 
numerator. In the next section, we will give an explicit construction of the BCJ box numer- 
ators for MHV amplitudes. For non-MHV amplitudes, however, another approach will be 
necessary, and we will return to this question in Section VIII. 
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IV. SELF-DUAL BOXES IN MHV AMPLITUDES 



In this section, we present a set of BCJ numerators for one-loop box diagrams in A/" = 4 
MHV amplitudes. There are two ingredients. The first is a conjecture relating one-loop 
amplitudes in self-dual gauge theory and in the MHV sector of the JV = 4 theory [45] . The 
second is the evaluation of the self-dual amplitudes using D-dimensional unitarity cuts [50] . 
Together, they allow for the construction of a linear map between BCJ box numerators and 
box-cut contributions to the self-dual amplitudes. 

Let us start with the remarkable conjecture of ref. [45]. Denoting by An-i the one-loop 
leading-colour partial amplitudes, it can be expressed as 

<5«(g)<^°(1^2^ . . . ,n+) = -2e(l - e)(47r)^ [<r(l,2, . . • ,n)]^^^^^ , (4.1) 

where D = 4 — 2e. On the left-hand-side, we have the one-loop amplitude in pure Yang- 
Mills theory with all-plus external helicities, which corresponds to the self-dual sector of the 
theory [51, 52]. The dimension shift D — )■ D + Aof the A/" = 4 amplitude is to be taken before 
integration, when the amplitude is expressed in terms of D-dimensional scalar integrals 
with m external particles, generically massive: 



The shift means that we substitute e -> e — 2 so that ^m^^- 

Although the expression above was conjectured to hold to all orders in e, we will only 
make use of the e — )■ limit. In this limit, we have D — > 8 as the dimension shift, and 
the right- hand- side of Eq. (4.1) is non- vanishing because the box integrals all diverge in 
the ultraviolet as (47r)^j|~^'^ ~ i/Qe. On the other hand, higher-degree scalar integrals 
(pentagons, hexagons, etc.) are ultraviolet finite in eight dimensions, and therefore will not 
contribute when e — > 0.^ Hence we find 



^'(Q)<r°(l+,2+,...,n+) = --^ <r'(l,2,...,n) , (4.3) 

boxes J4 



where the sum is over the kinematic coefficients of the scalar box integrals. 

We would now like to decompose the self-dual amplitude into an analogous sum over box 
contributions, so that we have the chance to identify the elements in that sum with the 
kinematic coefficients of the box integrals in the A/" = 4 amplitude. Such a decomposition 
of the self-dual amplitude is exactly what the unitarity method of ref. [50] provides. 

The self-dual amplitude, which is the one-loop amplitude in pure Yang-Mills theory with 
all-plus external helicities, has vanishing four-dimensional unitarity cuts. However, it can 
be constructed from D-dimensional box cuts, as shown in [50], 



^ n j— 3 j— 2 t— 1 

(4vr)^<r°(l+, 2^ • • • , n+) = -- 5^ 5^ E , (4-4) 



i=l j=i+l k=j+l l=k+l 



^ Even if pentagons, hexagons, etc., have a tensorial structure, the dimension shift still works if this structure 
does not lead to a pole in e. For instance, a pentagon integral whose numerator depends linearly on loop 
momentum does not contribute. 
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where we identify n + 1 ~ 1 and take e — )■ 0. Each box cut C is computed in the same manner 
as a four-dimensional maximal cut, except that the four-dimensional loop momentum is 
infinitely massive, i"^ = — £(_2^) = //^ — )■ oo, instead of being on-shell. In this limit of large 
11, the two cut solutions are given by 

£± = ±/i£ + 0(/i°), with = , 00^ = e^''P''K,,K2pK,^, (4.5) 

where we partitioned the external momenta into the four corners of the box cut C, so that 
Ki = pi+i + . . . + pj, and so on. Notice that, because of the Levi-Civita symbol in u^, if we 
permute the corners of the box the momentum i will be the same up to a change of sign, 
which will not affect the cut result. The contribution from a box cut C is 

Ci+i,... j|i+i,...,fe|fc+i,...,«|z+i,...,i = 2i lim A2 A3 Ai\ - , (4.6) 

where the A^. are the subamplitudes in each corner of the box. For instance, Ai is the 
subamplitude for the absorption of the external gluons i -|- 1 to j by a massive scalar of 
momentum For ^2, we have the absorption of the gluons j ' + 1 to by the same massive 
scalar, and so on; the mass of the scalar can always be taken as /x since the corrections from 
the gluons' momenta is 0{fi^). The leading behaviour of each Ar is linear in /i. We need 
these tree-level amplitudes up to six points for the cuts that we consider in this paper; they 
can be found in refs. [45, 46, 50, 53], and an all-multiplicity formula is given in ref. [54]. 

At this stage, we can already make an observation that will be crucial later on. Each 
box cut is symmetric under permutations of its corners, provided the order of the external 
gluons in each corner is preserved. That is, 

(4.7) 

The importance of this fact is that BCJ box numerators should also be symmetric under 
the permutation of its corners; otherwise, the Jacobi relations will require the presence of 
BCJ triangle numerators, which (we assume) are not required for A/" = 4 amplitudes. Let 
us also mention another useful property: 

Ci+l,...,i|i+l,...,A:|A:+l,...,;|«+l,...,i = (" l)"'"^*''*'^^'''''''^Qnv(i+l,...,j)|i+l,...,A:|A:+l,...,;|Z+l,...,i , (4.8) 

where inv(123 ■ ■ ■ m) is the inverse ordering m ■ ■ ■ 321, and # gives the number of external 
gluons in that corner (then # counts the three-point vertices which get inverted in Ar, while 
the other minus sign comes from the flip i — )■ —i). We will now give some examples, up 
to seven-point amplitudes. The six- and seven-point cases already illustrate all the issues 
involved at higher points. 

Four Points 

In the four-point case, we get 

(47r)2A^r°(l+,2+,3+,4+) 



V - ^MM (AO) 
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As is well known, this last expression does not depend on the ordering of the particles, so 
that Eq. (4.7) is satisfied. It follows from the dimension-shifting formula Eq. (4.3) that the 
box contributions to the D-dimensional A/" = 4 amplitude give 

A^=i - in'^n (1 9 A-\ fr rl|2|3|4,„ .1|3|4|2 rl|4|2|3\ . . . „x 



^4(1,2,3,4) = --(5«(g)Ci|2|3|4 = ^5'{Q)j^,^, (4.11) 



with 

[12] [34] 

2" ^-^^^i^i^i- ^^^(12)(34) 
where we defined 

Ci 23 4 = f'^^"'^^^ ^b2a2b3 Jb^a^bi Jbt^a^bi -|^2^ 

1|2|3|4 

and /4 is the box integral with the corresponding propagator structure. This well-known 
expression is already the full four-point MHV amplitude. To confirm that it is in BCJ form, 
we need to check that the Jacobi identities hold. The only type of Jacobi identity is the one 
connecting BCJ boxes to BCJ triangles: 

n4(l, 2, 3, 4) - n4(2, 1, 3, 4) = n;{[l, 2], 3, 4) . (4.13) 

In our case, 71,4(1, 2, 3, 4) is permutation symmetric. Therefore, the BCJ triangles vanish. It 
will still be true at higher points that the BCJ boxes are the lowest-order polygon required, 
as a consequence of the property (4.7). 



Five Points 

In the five-point case, we have 

(4vr)^A^f °(1 + , 2+, 3+, 4+, 5+) = -g (Cl2|3|4|5 + C'23|4|5|l + C'34|5|l|2 + C'45|l|2|3 + C'51|2|3|4) , 

(4.14) 

where 

ii|3|4|5 (I2)tr5(1234) ' ^ ' 

and similarly for the cyclic permutations. 

We see now how useful the decomposition (4.4) is; we can easily recognise to which box 
integral these different cuts correspond. The BCJ box contributions to the Af = 4 amplitude 
are then 



with 



BCJ box 



. 5 /n4([l,2],3,4,5) .(1+2)|3|4|5 , „ \ 

tg I C[i,2],3,4,5/4 +9 perms 1 , (4.16) 



n4([l,2],3,4,5) = -^5«(Q)s,2C,2|3|4|5 = '^'(Q)^^^^^^ ■ (4-17) 

The 10 permutations in Eq. (4.16) correspond to different choices of the two particles in the 
massive corner of the box. We also defined 

C[l 2] 3 4 5 = jb^bib^ Jb2a4b3 jb3,azb4 jb^aefi^ j^g-j 



10 



and 

r(l+2)|3|4|5 f 1 



This set of BCJ boxes coincides with the Carrasco- Johansson representation we discussed 
in Section III. 



Six Points 



The six-point case requires for the first time an additional ingredient. Consider^ 

i( 

(4.20) 



(47r)2A^f 0(1 + , 2+, 3+, 4+, 5+, 6+) = (2Ci23|4|5|6 + 2Ci2|34|5|6 + Ci2|3|45|6 + Cychc) 



where 

_ ^. (545(6|1+2|3] [51] [64]^ - g46(5|l+2|3] [61] [54]^) [56] 
^''I'l'l' " (12)(23)tr5(1456)tr5(3456) 

_ ^^. (5|1+2|6] (6|1+2|5] [12] [34] [56]^ 
t.i2|34|5|6- (12)(34)tr5(1256)tr5(3456) ' 

Cl2|3|45|6 = C'l2|45|3|6 • (4-21) 

There are three different cubic diagrams for the one-mass boxes with particles 1,2,3 in the 
massive corner, such as Ci23|4|5|6 • The numerators of those diagrams are 

n4([[l, 2], 3], 4, 5, 6), n4([[2, 3], 1], 4, 5, 6), n4([[3, 1], 2], 4, 5, 6) . (4.22) 

In terms of these numerators, we have 

1.8.n^^ ^4([[1,2],3],4,5,6) n4([[2, 3], 1],4, 5, 6) 

--0 (,C^jOi23|4|5|6 = , 

/ Si2'5l23 '523'5l23 

1 8 n4([[l,2],3],4,5,6) n4([[3,l],2],4,5,6) 

--d (^yj02l3|4|5|6 = \ , 

2 S12S123 S13S123 
^,S(n^r r24([[2,3], 1],4,5,6) n4([[3, 1],2],4,5,6) 

--0 (yjC231|4|5|6 = • l4.23j 

2 S23S123 S13S123 

Only two of these relations are independent, and one can verify the identity 

Cl23|4|5|6 + C'213|4|5|6 + C231|4|5|6 = , (4.24) 

using the explicit expressions in Eq. (4.21). Therefore, we need an additional relation in 
order to fix the three BCJ box numerators in terms of the box cuts. This relation is naturally 
the Jacobi identity: 

n4([[l, 2], 3], 4, 5, 6) + n4([[2, 3], 1], 4, 5, 6) + n4([[3, 1], 2], 4, 5, 6) = . (4.25) 
We may now solve for the numerators in terms of the cut boxes, and find 

n4([[l, 2], 3], 4, 5, 6) = -\8\Q)s^^ (s23Cl23|4|5|6 - Sl3C^213|4|5|6) • (4-26) 



The relative factors ensure that there is no over-count for the box Ci2|3|45|6 — C45I6I12 



2 

i 1 — 045|6|12|3- 
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The BCJ box numerators corresponding to the two-mass boxes can be obtained in a 
straightforward manner, 

n4([l,2],[3,4],5,6) = -\6\Q)s^2Sz^C,2m^^ ■ (4.27) 

Notice that, again due to the property (4.7) of the cut boxes, there is no difference between 
the so-called two-mass-hard (massive corners opposed) and two-mass-easy (massive corners 
adjacent) boxes, 

n4([l, 2], 3, [4, 5], 6) = n^{[l, 2], [4, 5], 3, 6) . (4.28) 

The BCJ box numerators obtained in this section coincide with the ones chosen by Yuan 
in ref. [47]; the present discussion may be regarded as a systematic derivation of these boxes. 
We will show below why it was not possible in ref. [47] to obtain the BCJ pentagon and 
hexagon numerators satisfying all the Jacobi identities: those numerators will depend on 
the loop momentum. Indeed, we will see that this dependence is fixed by the choice of box 
numerators. 



Seven Points 

At seven points, the simplest case is that of the BCJ numerators for three-mass boxes, 
n4([l,2], [3,4], [5, 6], 7) = -\6\Q)s^2S^^s^,C,2\^,\,^j . (4.29) 

For the two-mass boxes, one has to use the Jacobi identity as in the six-point case, with the 
result 

n4([[l,2],3], [4, 5], 6, 7) = -^5^(Q)si2S45 (s23Ci23|45|6|7 - ■531C312145I6I7) • (4.30) 

For the one-mass boxes, we have two different topologies of cubic diagrams, 

n4([[[l, 2], 3], 4], 5, 6, 7) and n4([[l, 2], [3, 4]], 5, 6, 7) . (4.31) 

If we consider all the permutions of the particles 1 to 4, these topologies give us 12 and 
3 different cubic diagrams, respectively. These 15 numerators satisfy 10 Jacobi identities 
among themselves, but only 9 of those are independent identities. On the other hand, after 
considering the property (4.8), there are 12 different corresponding box cuts, such as 

l.s.n^r n4([[[l,2],3],4],5,6,7) ^[[[3, 2], 1], 4], 5, 6, 7) 

--0 (^^/jOi234|5|6|7 = \ 

^ ■Si2Si23'5l234 S23S123S1234 

n4([[[2,3],4],l],5,6,7) n4([[[4, 3], 2], 1], 5, 6, 7) 

■S23S234S1234 "5345234-51234 

^ n4([[l,2],[3,4]],5,6,7) _ 
S12S34S1234 

These satisfy identities of the type 

Cl234|5|6|7 + C2134|5|6|7 + C'2314|5|6|7 + C'2341|5|6|7 = , (4.33) 
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which reduce the number of independent box cuts to 6. Together with the 9 Jacobi identities, 
the independent box cuts determine the 15 BCJ box numerators. The latter can be expressed 

as 

^4([[[1,2],3],4],5,6,7) = -^5^(Q)Si2^S34(s23Ci234|5|6|7 " ■S13C2134I5I6I7) 

+ (-513 + S23) (■5l4C'3214|5|6|7 ■S24C3124|5|6|7) + ■5i3Si4C23l4|5|6|7 ~ ■S23'524C'l324|5|6|7^ ; 

(4.34) 

and, obviously, 

n4([[l, 2], [3, 4]], 5, 6, 7) = n4([[[l, 2], 3], 4], 5, 6, 7) - n4([[[l, 2], 4], 3], 5, 6, 7) . (4.35) 



Higher Points 

We have seen in the previous examples that the box cuts C completely fix the BCJ box 
numerators. This will still be true for an arbitrary number of external gluons. The first 
step to see how this works is to notice that we can focus on one corner of the box cuts at a 
time. According to the expression (4.6), C factorizes into the four corners, which are only 
connected through the loop momentum cut solution. 

We want to identify the contribution from each corner of a box cut with a tree-level 
"BCJ current," consisting of external gluons and one off-shell leg connected to the loop 
momentum. Each "BCJ current" is a sum over diagrams with cubic vertices, whose kine- 
matic numerators satisfy Jacobi identities. The counting of independent BCJ numerators 
for rrir + 1 external legs is (m,. — 1)!. For instance, there were 6 such numerators in the 
one-mass box at seven points; recall the text below Eq. (4.31). 

To obtain the linear map between BCJ box numerators and box cuts C, we need a set of 
{rrir — 1)! independent box cuts. Such a set is the result of the Kleiss-Kuijf-type relations of 
the type: 

Cl23...m|a|/3|7 + C213...m\a\l3\^ + C231...m|a|/3|7 + . . . + C23...1m|a|/3|7 + C23...ml|o|/3|7 = , (4.36) 

where a, (3 and 7 are the sets of external gluons in the other three corners. These relations 
were proven in [37], where an explicit solution in terms of (m^, — 1)! independent box cuts 
was obtained, 

C<515'|a|^|7 = ( — 1)'''' ^ C'l(7|a|^|7 , (4.37) 
creOP{5T\j5') 

where the sum is over all permutations a that preserve the ordering of the elements in each 
subset 5'^ and 5'; 5'^ denotes the transposition of the elements in 5, and \5\ denotes the 
number of elements in 5. It was also shown in [37] that these relations are equivalent to the 
relations between self-dual partial amplitudes found in [55]. 

Let us mention the closely related work of ref. [56] , where certain BRST invariant kine- 
matic combinations found in superstring amplitudes have properties analogous to the box 
numerators constructed in this section. This may provide an alternative method to compute 
BCJ box numerators. 
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V. INTEGRAND REDUCTION 



In this section, we briefly discuss the problem of reducing scalar integrals algebraically 
at the integrand level. This will be important in order to compare the BCJ representations 
constructed in this paper with previously obtained representations of the amplitudes. 

Consider a scalar m-gon integral in D dimensions, given by 



which has massless propagators 



P^ = (i + p, + ...+p,^\ (5.2) 



and whose external momenta pi need not be null. With all external momenta in four 
dimensions and m > 5, we can partial fraction the integrand into a sum of scalar pentagon 
integrands 

1 ^ r((j5 + l) . . . + . . . J2|(j2 + 1) • • • J3|(j3 + 1) • • • J4|(j4 + 1) • • • Js) 

n"" p 2^ p p p p p ' ' 

where the coefficients r{. . .) depend on five ordered partitions of the set of external legs. By 
definition, r is cyclic in its arguments: r-(a;i|a2|«3|«4|a5) = ''"(a2|a3|«4|«5|«i)- Although the 
notation in Eq. (5.3) may seem needlessly lengthy, it will be crucial to be able to distinguish 
between different orderings of the external legs, so that, for example, r(12|3|4|5|6) and 
r(21|3|4|5|6) are not the same quantity. In particular, by r(21|3|4|5|6) we mean the coefficient 

of (•^^(^ + ^12)^(^ + ^123)^(^ + ^1234)^(^ + ^12345)^) ^ in the partial fractioning of the scalar 
hexagon integrand with external legs 213456, where pi2 = P1+P2 and so forth. We may also 
write something like r((l+2)3|4|5|6|7), which is the coefficient of the pentagon with corners 
(1+2+3, 4, 5, 6, 7) in the partial fractioning of a hexagon with corners (1+2, 3, 4, 5, 6, 7). 

This partial fractioning of scalar integrands was given explicitly in ref. [57] in terms of 
modified Cayley determinants. More recently, equivalent expressions were given in ref. [58] , 
which contains a more careful consideration of the dimensional regulator. 

In the present work, we prefer to use somewhat less cumbersome expressions in terms 
of traces of gamma matrices. As a simple example, consider the reduction of a hexagon 
integral, where we have 

r(12|3|4|5|6) = - ^^^(3456) 
^ I I I I ^ tr5(123456) ^ ' 

Eq. (5.4) can be proven straightforwardly by verifying the identity 

1 = r(12|3|4|5|6)Pi + r(l|23|4|5|6)P2 + r(l|2|34|5|6)P3 

+ r(l|2|3|45|6)P4 + r(l|2|3|4|56)P5 + r(61|2|3|4|5)P6 , (5.5) 

for example by considering the quantity tr5(12345£) and anti-commuting i around the trace. 
As a more complicated example, the reduction of an octagon integral requires the quantity 

n23l4l5fil7l8l = tr5(4(5+6)78) tr5(4(5+6)78) tr5((l+2+3)478) 

^ ''^ tr5(l(2+3)4(5+6)78)tr5((l+2)34(5+6)78)tr5((l+2+3)45678) ■ ^ ' 
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The three factors in this expression correspond to the three places one can insert an extra 
propagator in the pentagon specified by (123|4|56|7|8) while respecting the ordering of the 
external legs. Following these two examples, it is not difficult to guess the general expression, 
which follows by combining Eq. (5.4) with the results of ref. [57]: 

r(ai|«2|«3|«4|«5) = ll 11 (-—7- — -^1, (5.7) 

where {f3i, (32} = aj indicates all ways of splitting aj into two subsets f3i and that respect 
the ordering of aj, and the hat in the numerator indicates the absence of from the trace. 

It is important to recognise that the reduction of a scalar m-gon to scalar pentagons in 
Eq. (5.3) is an algebraic identity on the integrand, relying only on having four-dimensional 
external momenta and not on any integration identities such as vanishing of total derivatives. 
This will be useful in our discussion of the BCJ form of the six- and seven-point amplitudes, 
whose numerators depend on loop momentum; it will allow us to perform integrand reduction 
while leaving the numerators unspecified. 

It is also possible to reduce scalar pentagon integrals to scalar boxes, as in refs. [58-60]. 
Unlike reduction to pentagons, however, reduction to boxes is not an algebraic identity; it 
throws away parity-odd integrands that integrate to zero and is only valid up to terms of 
0(e). Nevertheless, we will find it useful to compare with known expressions for amplitudes, 
given as linear combinations of scalar box integrals. The reduction from scalar pentagons 
to scalar boxes is then achieved using the reduction coefficient 

.1oloul.^ _ tr(34512345) -4(pi ■P2)p^p1p^ 
n[iZ\6m^) - 4p2p2p2p2p2_t,(i234512345) " ^^"^^ 

Finally, one can explicitly reduce m-gons to boxes by first reducing to pentagons, yielding 

4 

r{ai\a2\a3\a^) = Y, ni- ■ ■ If^M ■ ■ M- ■ ■ \W2\ ■ ■ ■) ■ (5.9) 

i=i {/3i,/32}=«j 



VI. THE SIX-POINT MHV AMPLITUDE 



We are now in a position to compute the numerators of the six-point one-loop MHV 
amplitude. In fact, the computation is remarkably straightforward if we follow the guidance 
from the five-point case. We begin with the box numerators discussed in Section IV, namely 

^4([[1, 2], 3], 4, 5, 6) = -^6^{Q)Si2 (s23Cl23|4|5|6 - Si3C2i3|4|5|6) , (6.1) 
n4([l,2],[3,4],5,6) = -^5«(g)Si2S340i2|34|5|6 • (6.2) 

These box numerators do not depend on loop momentum. More generally, we do expect the 
numerators to depend on loop momentum; indeed, it is reasonable to guess that the m-gon 
numerator will have up to m — 4 powers of loop momentum. 

Our first task is to determine the pentagon numerators n^{a,b,c,d,e; i).^ Since we are 
dealing with a six-point amplitude, one of the legs of the pentagon is massive; we will treat 

^ Recall that the subscripted 5 refers to the pentagon, regardless of what multiplicity amplitude we are 
discussing. 
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both massive and massless legs on the same footing. We begin by following the integrand 
oxidation procedure we described in Section III - namely, we move one leg around the loop. 
In detail, we consider these equations: 



715(0, b, c, d, e; i) — n^{b, a, c, e; t) = n4^([a, b], c, d, e) , (6.3) 
n^i^b, a, c, d, e; i) — n^i^b, c, a, d, e; £) = riiib, [a, c], (i, e) , (6.4) 

n^{b, c, d, a, e; i) — n^{b, c, d, e, a; £) = "0.4(6, c, [a, e]) . (6.5) 

Notice that we have moved the leg a around the loop, so we expect that n^i^b, c, d, e, a; £) is 
related to n^i^a, b, c, d, e; £). Indeed, the cyclic permutation property of the pentagon is 

71,5(6, c, d, e, a; i) = n^{a, b, c, d,e]i — pa) ■ (6.6) 

Therefore, by summing our system of equations, we learn that 

n^{a, b, c, d, e; i) — n^{a, b, c, d,e;i — Pa) = x^{a, b, c, d, e) , (6.7) 

where 

^5(0, b, c, d, e) = n4([a, b], c, d, e) +n4(6, [a, c], rf, e]) +n4(6, c, [a, d], e) +n4(6, c, (i, [a, e]) . (6.8) 

Now, let us consider expanding 71,5(0, b, c, d, e; £) in the loop momentum i: 

n^{a, b, c, d, e; £) = ^0,5(0, b, c, e) + ^11,5(0, b,c,d,e) ■ i + i ■ n2,5(a, 6, c, ci, e) ■ £ + ■ ■ ■ (6.9) 

The structure of Eq. (6.7) is basically a finite difference equation. Since 0:5(0, 6, c, d, e) is 
independent of loop momentum (after all, it is a sum of boxes), we see that 115 can depend 
on at most one power of loop momentum (n2,5 = 0, etc.), consistent with our expectation. 
We now find that 

Pa ■ ni^5{a, b, c, d, e) = ^5(0, b, c, d, e) . (6.10) 

This is not all we know about the numerator ni 5. From Eq. (6.3), whose right-hand side 
is independent of the loop momentum, and from the vanishing of 722,5, we conclude that 
ni,5 is permutation symmetric. This is enough to determine ni 5 by cyclically permuting 
Eq. (6.10): 

Pb ■ ■'^1,5(0-, b, c, d, e) = X5{b, c, d, e, a) , (6.11) 
Pc ■ ni^^^a, b, c, d, e) = ^5(0, d, e, a, b) , (6.12) 
Pd ■ ^1,5(0, b, c, d, e) = X5{d, e, a, b, c) . (6.13) 

We can therefore reconstruct the numerator as 

ni^{a, b, c, d, e) = ^5(0, b, c, d, e) + ^5(6, c, d, e, a) 

+ cu^ 0:5(0, d, e, a, 6) + X5(ci, e, a, 6, c), (6.14) 



where 



e(-,6,c,d) e(a, -,0,^) 

= ^ — ; 3V ' '^fe = — I ' ^tc. (6.15) 

e(a,o, c, a) e[a,b,c,d) 
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The vectors {ua,0Jb,(^c,^d} form a four-dimensional basis dual to {pa,Pb,Pc,Pd}- 

To determine the part of the pentagon numerator that is independent of loop momentum, 
we use Jacobi moves to rearrange the legs from a starting position to its reverse. Then we 
can use the mirror symmetry of the pentagon to deduce a simple equation. Prom the Jacobi 
moves, we deduce that 



n5(a, b, c, d, e; £) — n5(c, b, a, e, d; £) = y^{a, b, c, d, e) , 



(6.16) 



where 



?/5(a, 6, c, d, e) = n4([a, 6], c, rf, e)-\-ni{b, [a, c], d, e)-|-n4([6, c], a, d, e)-|-n4(c, 6, a, [d, e]) . (6.17) 
The symmetries of the pentagon require that 




•H- -^5(0,6, c, d,e;-£ + pd+Pe), (6.18) 



allowing us to determine 710,5 in terms of the boxes we constructed from the self-dual theory 
and the linear pentagon produced with the oxidation procedure, 

no,5(a,&,c,(i,e) = ^ (1/5(0, 6, c, ci, e) - {pd+Pe) ■ ni^5{a,b,c,d,e)) . (6.19) 

We have now fully constructed the pentagon numerator. It is straightforward to construct 
the hexagon numerator using the same methods. First, we move one leg around the loop 
and use cyclicity to discover that 



naia, b, c, d, e, /; £) - 715(0, b, c, d, eJ-J- Pa) = ^6(0, 6, c, d, e; £) , 



(6.20) 



where 



XQ{a,b,c,d,e]£) = n5{[a,b],c,d,e, f;£) +n5(6, [a,c],d,e, f;£) +715(6,0, [a,d],e, f;£) 

+ n^(b, c, d, [a, e], /; £) + n^{b, c, d, e, [a, /];£). (6.21) 

Notice that ,x'6. which is a sum of pentagons, may now depend on the loop momentum. Since 
there is at most one power of £ in Xg, the hexagon numerator can have at most two powers 
of £. However, a remarkable cancellation takes place. The linear combination of pentagons 
in xq is such that it does not depend on the loop momentum after all, a fact which boils 
down to our choice of box numerators. Therefore, the hexagon numerator depends on only 
one power of loop momentum. 



7i6(a, 6, c, d, e, /; £) = nofi{a, b, c, d, e, /) + ni^e{a, b, c, d,eJ)-£. 



(6.22) 
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Using this expansion in Eq. (6.20), we find an equivalent expression involving ni^ alone: 

Pa ■ ^i,6(a, b, c, d, e, f) = X(i{a, b, c, d, e, f) . (6.23) 

Using cyclicity of ni g it is straightforward to find a solution: 

n^fi{a, b, c, d, e, /) = xe{a, 6, c, d, e, /) + x^ib, c, d, e, /, a) 

+ xe{c, d, e, /, a, h) + uj^ xe{d, e, /, a, b, c) . (6.24) 

Another non-trivial cancellation occurs when we perform an integrand reduction to pen- 
tagons, if we now think about a partial amplitude. The linear loop momentum dependence 
of the hexagon numerator cancels against that of the pentagon numerators: 

f h A f\ ( u \A\ \f\ ^ ^i,5([a,6],c,c?,e,/) 
ni,6(a, 0, c, d, e, /) r{ab\c\d\e\f) H = . (6.25) 

^ab 

The consistency between the expressions (6.24) and (6.25) shows the tight connection be- 
tween the system of Jacobi identities and the structure of integrand reductions. 

It remains to determine no,6- In the case of the pentagon, we used inversion symmetry 
to find the analogous part of the numerator. But for the hexagon this is no longer of any 
help. The reason is simply that riQ^a, b, c, d, e, /; i) = +nQ{f, e, d, c, b, a; —i) since we are now 
reversing the orientation of six antisymmetric cubic vertices. Thus, the inversion symmetry 
leads to an equation of the form 

2ni,6(«, b, c, d,e,f)-i = y^ia, b, c, d, e, /; i) , (6.26) 

where i/q is a linear combination of pentagons. In fact, the procedure we have been following 
of systematically exploiting the Jacobi relations and the symmetries of the loop diagrams 
cannot fix the scalar part of a hexagon numerator. To see this, suppose we have the correct 
scalar hexagon numerator no,6 in our hands. Consider adding a function S{a, b, c, d, e, /) of 
kinematic variables to this numerator: 

hofi{a, b, c, d, e, /) = no,6(a, b, c, d, e, /) + S{a, b, c, d, e, /) . (6.27) 

The Jacobi relations, and the cyclicity and inversion symmetries of the hexagon, are un- 
changed provided that S" is a permutation-symmetric function. To fix this ambiguity we 
must add more physical information. We choose to do so simply by requiring that our 
numerators reproduce the known MHV amplitude [43, 62]. We use integrand reduction 
to boxes; in particular, it is known that the two-mass-hard box coefficient vanishes. The 
reduction is straightforward, with the result that 

no,6(a, b, c, d, e, /) r{ab\cd\e\f) + ^' ^' ^' ■^K {a+b\cd\e\f) 

Sab 

+ ^' riab\c+d\e\f) + d],ej) ^^_ ^^^8) 

Scd Sabred 

This uniquely fixes the scalar part of the hexagon numerator. The Jacobi relation to the 
scalar pentagon is satisfied. 

In this way, we have built a set of numerators which satisfies all the Jacobi relations. In 
the construction, we used the absence of the two-mass-hard boxes. To verify our procedure, 
we numerically checked that the two-mass-easy and one-mass box coefficients take their 
known values. 



18 



VII. THE SEVEN-POINT MHV AMPLITUDE 



Essentially the same techniques suffice to compute the seven-point, one-loop, MHV am- 
plitude. A new element to emerge in this case is that there is a family of numerators for 
the heptagon with the property that the associated boxes are those computed from the self- 
dual boxes. In other words, there are generalised gauge transformations among colour-dual 
numerators at seven points that preserve the box numerators. 

The computation starts with the boxes. There are three kinds of box numerators relevant 
in this computation; namely, the one-, two- and three-mass boxes. As we discussed in 
Section IV, we may take the numerators of these boxes to be 

n4([l,2],[3,4],[5,6],7) = -^5'(g)Si2S34S56Cl2|34|56|7 , (7.1) 
n4([[l, 2], 3], [4, 5], 6, 7) = -^5®(Q)si2S45 (s23Cl23|45|6|7 - S31C312I45I6I7) , (7.2) 

for the three- and the two-mass boxes, respectively. All one-mass boxes can be computed 
from 

'^4([[[1,2],3],4],5,6,7) = -^5^((5)S12 1^534 (s23Ci234|5|6|7 " S13C2134I5I6I7) 

+ (513 + S23) (si4C32i4|5|6|7 -524^3124151617) + ■5l3'5l4C*2314|5|6|7 ~ ■S23'524Cl324|5|6|' 

(7.3) 

Given these boxes, we can construct pentagon numerators in complete analogy with the 
hexagon case. Thus, the pentagon numerator, 

n^{a, b, c, d, e; i) = no^^^a, b, c, d, e) + ^1^5(0, b, c,d,e) ■ i , (7.4) 

is given by 

ni^{a, b, c, d, e) = 0:5(0, b, c, d, e) + ^5(6, c, d, e, a) 

+ co^ X5(c, d, e, a, b) + X5{d, e, a, b, c) , (7.5) 
2no,5(a, b, c, d, e) = 1/5(0, b, c, d, e) - pde ■ ^1,5(0, b, c, d, e) , (7.6) 

where 

x^^a, b, c, d, e) = n^da, b], c, d, e) + n^^ib, [a, c], d, e) + n^ip, c, [a, rf], e) + ^4(6, c, rf, [a, e]) , 

(7.7) 

7/5(0, 6, c, e) = n4([a, 6], c, d, e) + n4(6, [a, c], (i, e) + n^{\b, c], a, e) + ^4(0, 6, a, [rf, e]) , 

(7.8) 

and, just as in Section VI, we define 

e(a,o, c, dj e(a, 0, c, aj 

In the case of the six-point amplitude, the hexagon numerator contained at most one 
power of the loop momentum. This is no longer the case at seven points. By walking leg a 
around the hexagon, we find that the hexagon numerator 

n6(a, 6, c, d, e, /; t) = nofi{a, b, c, d, e, /) +nifi{a, b, c, d, e, f)-i+i-n2,6{a, b, c, d, e, f)-i (7.10) 



19 



satisfies 



2PaM ni^fiia, b, c, d, e, /) = x^Q^a, b, c, d, e, /) , (7.11) 
Pa ■ nifi{a, b, c, d, e, /) = Xofi{a, b, c, d, e, f)+Pa- ?^2,6(a, b, c, d, e, f) ■ Pa , (7.12) 

where x^^a, b, c, ci, e, /; i) = xofi{a, b, c, d, e, /) + xi^^^a, b, c, c?, e, /) ■ £ is a sum of pentagons 
obtained in tlie integrand oxidation algoritfim. In detail, 

xeia, b, c, d, e, /; £) = n^^a, b], c, d, e, /; i) + n^{b, [a, c], rf, e, /; £) 

+ 725(6, c, [a, ci], e, /; £) + ^5(6, c, rf, [a, e], /; ^) + rig (6, c, d, e, [a, /];£). (7.13) 

Tfie key difference relative to tlie six-point construction is tliat, in contrast to tlie discussion 
below Eq. (6.21), in this case the quantity xi^^a, b, c, d, e, /) does not vanish. 

It is straightforward to solve for n2fi{a, b, c, d, e, /). Notice that, without loss of generality, 
we can take n2,6 to be a symmetric tensor. Moreover, since the pentagon numerator is at most 
linear in i, all the Jacobi relations involving n2,6 must be trivial. Therefore, n2fi{a, b, c, d, e, /) 
is a permutation symmetric function of the legs a, b, c, d, e and /. A symmetric tensor is 
determined if we know its projections onto a basis of four vectors, which we take to be pa, 
Pb, Pc and pd- (Note that these momenta may be off-shell, corresponding to a massive leg.) 
Thus, we need to construct a tensor satisfying 



2pa • n2fi{a, b, c, d, e, /) = xi,6(a, b, c, d, e, /) , (7.14) 

2ph ■ n2fi{a, b, c, d, e, /) = Xifi{b, c, d, e, /, a) , (7.15) 

2pc ■ n2fi{a, b, c, d, e, /) = xi^q{c, d, e, /, a, b) , (7.16) 

2pd ■ n2,6(a, 6, c, rf, e, /) = Xi,6(ci, e, /, a, 6, c) . (7.17) 

It is non-trivial that these equations are consistent with the symmetry of ^25. Using 
Eq. (7.14) and Eq. (7.15), we see that 

2pa ■ ^2,5(0, &, c, d, e, f) ■Pb = Pb- xifi{a, b, c, rf, e, /) , (7.18) 

2pfe ■ n2fi{b, c, ci, e, /, a) ■ pa = Pa ■ xifi{b, c, e, /, a) . (7.19) 

Therefore consistency requires 

Pb ■ xifi{a, b, c, d, e, f)=Pa- xifi{b, c, c?, e, /, a) . (7.20) 



We have verified this equation using the explicit self-dual boxes. In view of this consistency, 
the solution for the numerator 77-2 g is simply 

"-2,6 6, c, d, e, /) = ul Pa ■ xi,6(a, 6, c, c?, e, /) + ■ xi,6(&, c, d, e, /, a) 

+ uj^ uf^ Pc ■ xi,6(c, d, e, /, a, h) + ■ Xx,%{d, e, /, a, 6, c) 

+ (^a^6 + ^h^l) Pb ■ xi^eia, b, c, d, e, /) + (w^^ + w^^;;) ■ Xi^^ia, b, c, c?, e, /) 

+ i^a^d + ^d^a) Pd ■ xi^&{a, b, c, d, 6, /) + (w^w^ + w^w^ ) " a;i,6(&, c, rf, e, /, a) 

+ i^b^d + ^d^b) Pd ■ xi,6{b, c, d, e, /, a) + (w^c^^ + uj^^u'^) pd ■ xi,q{c, d, e, /, a, 6) . 

(7.21) 
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With ^2,6 in hand, we now tTirn our attention to ui^q, the part of the hexagon numerator 
that is hnear in loop momentum. It is determined by 

Pa ■ "'i,6(a, b, c, d, e, /) = xo,6(a, b, c, d, e, f)+Pa- ^2,6(0, 6, c, d, e, /) • p„ 

= Xo,6(a, 6, c, e, /) + ^a;i,6(a, &, c, c^, e, /) • Pa . (7.22) 

We have met equations hke this several times already; in this case, there is a slight novelty 
since ni^Q{a, b, c, d, e, f) is not cyclically symmetric. Instead, the linear numerator satisfies 

ni,6(6, c, d, e, /, a) = nifi{a, b, c, d, e, f) - 2pa ■ n2fi{a, b, c, d, e, /) 

= ni,6(a, 6, c, ci, e, /) - a;i,6(a, 6, c, d, e, /) . (7.23) 

Taking this into account, the unique solution is 

riifiia, b, c, d, e, /) ^u^ria + rjb + uj^ric + m , (7.24) 

where 

Va = xofiia, b, c, d, e, f)+Pa- (^Xifi{a, b, c, d, e, /)^ , (7.25) 
Vb = Xo,6{b, c, d, e, /, a) + Ph ■ (^Xifi{b, c, d, e, /, a) + Xi^eia, b, c, d, e, /)^ , (7.26) 
?7c = a;o,6(c, (i, e, /, a, b) + Pc ■ (^xifi{c, d, e, /, a, b) + xifi{b, c, d, e, /, a) 

+ xi,6(a, 6, c, d, e, /)^ , (7.27) 
% = xo,6{d, e, /, a, 6, c)+pd- (^xifi{d, e, /, a, b, c) + xi,6(c, d, e, /, a, &) 

+ xi^e{b, c, 0?, e, /, a) + Xi^e{a, b, c, d, e, /)^ . (7.28) 

We have now constructed all of the hexagon numerators other than the scalar hexagon, 
riofi. Our methods do not fix this part of the hexagon numerator at this stage, so we 
will continue to the heptagon. As usual, the oxidation algorithm yields a simple equation 
satisfied by the numerator 

717(0, 6, c, d, e, /, g; £) - njib, c, d, e, /, g, a; £) = 2:7(0, 6, c, d, e, /, g; £) , (7.29) 

where ^7(0, b, c, d, e, /, g; i) is given by the familiar sum over hexagons: 

X7{a,b,c,d,eJ,g;i) = ne{[a,b],c,d,ej,g;i)+ne{b, [a,c],d,ej,g;i) 

+ 716(6, c, [a, d] , e, /, g; i) + n&{b, c, d, [a, e] , /, g\ £) 
+ neib, c, d, e, [a, /] ,g;i)+ nQ{b, c, d, e, /, [a, g^] ; £) . (7.30) 

Since uq contains up to two powers of loop momentum, in principle rij could contain up to 
three powers. However, we find that this does not occur since the terms in xj which contain 
two powers of loop momentum cancel out. We therefore expand 

7i7(a, 6, c, d, e, /, g) = 710,7(0, 6, c, d, e, /, g) + 711,7(0, 6, c, d, e, /, g') • ^ 

+ e-n2,7{a,b,c,d,ej,g) (7.31) 
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and find that the various different parts of the numerator satisfy 

2pa ■ n2j{a, b, c, d, e, /, ^f) = Xij{a, 6, c, d, e, /, ^f) , (7.32) 

Pa ■ nij{a, b, c, d, e, f, g) = xoj{a, b, c, d, e, f, g) + ]^xij{a, b, c, d, e, f, g) ■ Pa ■ (7.33) 

These are exactly the equations the hexagon numerators satisfy. They may be solved by 
precisely the same methods as in Section VI; although n2j{a, b, c, d, e, /, g) is not permutation 
symmetric, it is still cyclically symmetric, which is sufficient. Meanwhile, the cyclic shift of 
riij is in complete analogy with Eq. (7.23), namely 

nij{b, c, d, e, /, g, a) = nij{a, b, c, d, e, /, ^f) - Xij(a, 6, c, d, e, /, ^f) . (7.34) 

At this point, we have not constructed the part of the hexagon numerator which is indepen- 
dent of loop momentum, so we do not have xqj at our disposal. We do have xij and may 
therefore solve for n2j. It is 

n^7(a, 6, c, d, e, /, g) = Pa ■ xij{a, b, c, d, e, /, g) + pb ■ Xi^jib, c, d, e, /, g, a) 

+ Pc ■ a;i,7(c, d, e, /, g, a, b) + uj""^ pd ■ xij{d, e, /, g, a, b, c) 

+ i^a^b + ^b^a) Pb ■ xijia, b, c, d, e, /, g) + {uj^u'^ + u^u^) Pc ■ Xij{a, b, c, d, e, /, g) 

+ i^a^d + ^d^a) Pd ■ xij{a, b, c, d, e, /, g) + (w^w^ + uj^ul) Pc ■ xij{b, c, d, e, /, g, a) 

+ i^b^d + ^d^b) Pd ■ xi,7{b, c, d, e, /, g, a) + {u'^u^ + uJd^c) Pd ■ xi,7{c, d, e, /, g, a, b) . 

(7.35) 

A non-trivial property satisfied by this numerator is that upon integrand reduction to pen- 
tagons, it cancels against n2,6; this property is analogous to what we saw at six points, in 
Eq. (6.25). In fact, if we impose that on integrand reduction to pentagons, the terms linear 
in i also cancel, we can determine riij. It is given by 

. , , . ^ -1 f nifi{[a,b],c,d,ej,g) , , , , x 
nij{a,b,c,d,e, f,g) = . , , , , r{a+b\cd\e\f\g) 



r{ab\cd\e\f\g) \ Sab 
+ ^^-^^"'^'[^'^]'^'^'^V(a5|c+^|e|/|,)) + ^i.5(K&]Jc,ci],e,/,,)\ ^^^^^ 

Scd Sab-S^d J 

It is non-trivial that this expression for nij{a,b,c,d,e, f, g) has the property that all the 
Jacobi identities are satisfied. 

Encouraged by this success, we recall that Eq. (7.33) shows that ni 7 is determined by 
the scalar parts of the hexagons; turning this around, the scalar hexagons are constrained 
by knowledge of riij. Let us begin by understanding how much freedom there is at this 
point in the scalar hexagon numerators. Since we know all the scalar pentagons, we can use 
Jacobi moves to bring a given scalar hexagon to a canonical order. We choose to place the 
unique massive leg of the hexagon in the first entry of the function; we may order the rest 
of the legs numerically. There are 21 choices of the massive leg, so as it stands there are 
21 functions of the form nofi{[a,b],c,d,e, f, g) to be determined; all other scalar hexagons 
can be determined from these 21 functions and the scalar pentagons. The structure of the 
equations relating the scalar hexagon numerators no,6 to the linear heptagon numerator riij 
is as follows. Define 

si(a, 6, c, d, e, /, g) = Pb ■ ?^i,6([a, c], d, e, /, g, b) + p^c ■ ^i,6([a, d],e, /, g, 6, c) 

+ Pbcd ■ ni^fi{[a,e]J,g,b,c, d) + Phcde ■ ni^(i{[aj],g,b,c, d, e) 

+ Pbcdef nifi{[a, g],b, c, d, e, f) (7.37) 
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and 

S2{a,b,c,d,eJ,g) = Pb ■ n2,6i[a, c], d, ej, g,b) -pb + Pbc ■ n2,6{[a,d],ej, g,b,c) ■ pbc 

+ Pbcd ■ n2fi{[a,e],f,g,b, c, d) -Pbcd+Pbcde ■ n2fi{[a, f], g,b,c,d,e) ■ Pbcde 
+ Pbcdef ■ n2fi{[a, g],b, c, d, e, f) ■ pbcdef ■ (7.38) 

These quantities simply account for cycling the hexagon numerators to a canonical order. 
We then find 

nofiila, 6], c, d, e, f, g) + no,6([a, c], d, e, /, g, b) + no,6([a, rf], e, /, g, b, c) 
+ nofi{[a, e], f, g, b, c, d) + no,6([a, /], ^, b, c, d, e) + no,6([a, g\,b, c, d, e, /) 

= Pa ■ ^1,7(0, c, d, e, /, ^) - -Pa ■ xij{a, b, c, d, e, f, g) - Si(a, 6, c, e, /, «/) 

- S2ia,b,c,d,eJ,g) . (7.39) 

Note that the unknowns appearing in this equation include all possible choices of scalar 
hexagon numerator whose massive leg includes particle a. There can therefore be at most 
seven independent equations, corresponding to the different possible choices for leg a. In 
fact, only six of these equations are independent. We may therefore use these equations to 
solve for the six different numerators whose massive leg include a given particle, for example, 
particle 7. Thus, 

no,6([a,7],b, c, d,e,f) 

= Pa ■ nij{a, b, c, d, e, /, 7) - -pa ■ xij{a, b, c, d, e, /, 7) - si(a, b, c, d, e, /, 7) 

- S2{a,b,c,d,e, f, 7) - nofi{[a,b], c, d,e,f, 7) - no,ei[a, c], d, e, f, 7, b) 

- no,e{[a, d],e, /, 7, 6, c) - no,6([a, e], /, 7, 6, c, d) - no,6([a, /], 7, 6, c, e) . 

(7.40) 

There are 15 scalar hexagon numerators still to be fixed. We will fix them by integrand 
reduction. But before we do so, we must compute the family of scalar heptagon numerators 
which correspond to the space of hexagon numerators we are considering. 

Given a scalar hexagon numerator, it is straightforward to compute the scalar heptagon 
numerator. We use the refiection symmetry of the heptagon to do so; one convenient method 
is to note that 

717(0, 6, c, d, e, /, g; £) - 717(0?, c, 6, a, g, /, e; £) = t/7(a, 6, c, d, e, /, g; i) , (7.41) 

where yr is the following sum over hexagon numerators 

yr{a,b,c,d,ej,g;i) = nG{[a,b], c, d, ej, g; i) + n(i{b, [a,c],d,e, f,g;i) 

+ n6ib,c, [a,d],ej,g;i) + n6{[b,c],d,a,e, f,g;i) + rieic, [b,d],a,e, f,g;i) 
+ nedc, d], b, a, e, f, g; i) + n6{d, c, b, a, [e, f],g; i) + n6{d, c, b, a, /, [e, g]-J) 
+ nQ{d,c,b,a,[f,g],e-e). (7.42) 

The symmetries of the heptagon require 

7i7(d, c, 6, a, g, /, e; £) = -717(0, 6, c, d, e, /, g-; -i + Pe+Pf + Pg) , (7.43) 
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so that we can deduce a very simple expression for the scalar heptagon numerator noj in 
terms of the scalar hexagons and other known quantities: 

noj{a,b,c, d, e, f,g) 

= ^ {yoj{a, b, c, d, e, f, g) - Pefg ■ nij{a, b, c, d, e, f, g) - Pe/g ■ n2j{a, b, c, d, e, f, g) ■ Pefg) ■ 

(7.44) 



So we have constructed a 15-parameter family of numerators for the seven-point amplitude 
that satisfies all the Jacobi relations. We now want to exploit the freedom in our construction 
to make contact with the known expressions for the MHV boxes on integrand reduction. 
One simple property of the box coefficients is that the three-mass-hard boxes all vanish 
[43, 62]. We find that this requirement leads to 13 independent equations of the form 



= nojia, b, c, d, e, /, g) r{ab\cd\ef\g) + ^o.6([a, b],c, d, e, f, 9) 

Sab 

+ ^°'^^"'^'[^'^]'^'^'^V (a&|c+^|e/|,) + ^^MM^Al^^r(a6|crf|e+/|^ 



Scd 



+ ^°-^^["'^]'[^'^]'^'^'^V (a+6|c+c^|e/|,) + ^' ^' ^V (a+6|c^|e+/|,) 



+ 



SabScd 

wo,5(a,6, [c, d], [ej],g) 

ScdSef 



r{ab\c+d\e+f\g) 



SabSef 

n4([a,6],[c,(i],[e,/],^) 

SabScdSef 



(7.45) 



This is a linear system of equations which may be solved to fix 13 of the remaining unknown 
scalar hexagon numerators. Two numerators remain free. Indeed, we have more physical 
information at our disposal: we know that the two-mass-hard boxes vanish, while the two- 
mass-easy and one-mass boxes take on particular values. However, all of these requirements 
are automatically satisfied. We are therefore left with two functional parameters which may 
be chosen arbitrarily. This completes our construction of the seven-point amplitude. 



VIII. BEYOND MHV AMPLITUDES 



In the previous sections of this paper, we have exploited information originating in the 
self-dual sector of Yang-Mills theory as the basis for building complete sets of colour-dual 
numerators. Since the dimension- shifting formula is limited (as far as we know) to the 
MHV amplitudes, we present here a more general construction. We determine a linear map 
that fixes the scalar part no of the colour-dual numerators from unitarity cuts [43, 62], 
subject to the assumption that all i dependence in our numerators cancels upon reduction 
to pentagons, as happened in the MHV case. The fact that this linear map has solutions 
is due to nontrivial identities between the unitarity cuts. The i dependence is then fixed 
by the scalar part using the integrand oxidation methods described in previous sections. In 
the five-point case, the linear map is simple enough to be solved analytically; for six and 
seven points, we outline its construction and find solutions. We see no obstruction to the 
extension of this method to higher multiplicities. 
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Five Points 



At five points, there are only one-mass boxes available. In the BDDK presentation of the 
amplitude (i.e., after integrand reduction to boxes, throwing away the parity-odd terms) the 
box coefficients are given by 

vim; 4 ^^^(12) (23) (34) (45) (51) ^ ' 

Our goal is to construct a set of colour-dual numerators which integrand-reduce to these 
one-mass boxes. In other words, we wish to express the integrand of the colour-ordered 
amplitude as 





'^n4^{[a,b],c,d,e) ei , (8.2) 



where the sum is over the five cyclic permutations of (1,2,3,4,5). We will refer to this 
presentation of the integrand as the colour-dual presentation. Notice that we have assumed 
that the pentagon numerator is independent of loop momentum; this is consistent with our 
general expectations but can be thought of as an Ansatz for the purpose of the present 
discussion. 

Now, the colour-dual presentation must equal the BDDK presentation on integrand re- 
duction, which leads to 

,/ , , , N / N / , , , N '^sfl, 2, 3, 4, 5) — 77.5(2, 1, 3, 4, 5) 
12 3 4 5 =^5 1,2,3,4,5 r 12 3 4 5 + ' ' ' ' ^ 5^^^^^ _ g_3 

S12 

This equation is not especially transparent, but we may derive a simpler equation as follows: 
switch legs 1 and 2 in Eq. (8.3) to find 

- t 6(21|3|4|5) = ^5(2, 1, 3, 4, 5) r(21|3|4|5) + »5(2, 1, 3, 4, 5) - ^5(1, 2, 3, 4, 5) ^ 

and add this to Eq. (8.3). The result is obviously 
n5(l,2,3,4,5)r(12|3|4|5) + ^5(2, 1,3,4,5) r(21|3|4|5) 



^(5«(Q)tr((l+2)345) (^^2) (23) (34) (45) (51) + (21) (13) (34) (45) (52)) 
= X({1,2},3,4,5). (8.5) 

This equation is a kind of analogue of the Jacobi equation; the two terms on the left differ 
by an adjacent transposition. We may therefore consider moving one leg, for example 
leg 1, around the loop. In this way, we will obtain a simple equation for the numerator 
715(1,2,3,4,5). It is straightforward to handle the appearance of the reduction coefficient 
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factors. We find 

r(21|3|4|5)r(31|4|5|2)r(41|5|2|3) 



715(1,2,3,4,5) r(12|3|4|5) -r(51|2|3|4) 



r(13|4|5|2) r(14|5|2|3) r(15|2|3|4) 



Y(U 9\'^A^\ Y(U -XX A ^ o /(21|3|4|5) r(21|3|4|5)r(31|4|5|2) 
= 2}, 3, 4, 5) - 3}, 4, 5, 2) + ^({1, 4}, 5, 2, 3) ,(14|5|2|3) 

Yr/I 9 3 , /(21|3|4|5)r(31|4|5|2)r(41|5|2|3) 

A 5|, 2, 3, 4j^^^3|4|^|2) r(14|5|2|3) r(15|2|3|4) ' ^^"^^ 
We have checked that the solution of this equation is the famihar Carrasco- Johansson nu- 

Thus, we have recovered the known colour-dual set of numerators without starting from any 
particularly nice set of boxes. 



Six Points 



This method generalises to higher-point amplitudes. At six points, we can again match 
the colour-dual presentation of the amplitude to all of the box cuts: 

2-mass hard: - i 6(12|34|5|6) = no,6(l, 2, 3, 4, 5, 6) r(12|34|5|6) 

+ ^°'^^[^'']''''''''V (1+2|34|5|6) + ^o,5(l,2,[3,4],5,6) ^^^2^3^^^^^^^ 



+ 



S12 

no,4([l,2],[3,4],5,6) 
S12-S34 



S34 



(8.8) 



2-mass easy: - i 6(12|3|45|6) = no,6(l, 2, 3, 4, 5, 6) r(12|3|45|6) 



+ !5«(M.MlM,(i+2|3|45|6) + !^m(1iMMLS,(i2|3|4+5|6) 



+ 



SI2 

no,4([l,2],3,[4,5],6) 

S12S45 



S45 



(8.9) 



1-mass: - i 6(123|4|5|6) = no,6(l, 2, 3, 4, 5, 6) r(123|4|5|6) 



+ "°-^^t^''^'''''''' r((l+2)3|4|5|6) + 'V (l(2+3)|4|5|6) 

^ no,4([l,[2,3]],4,5,6) ^ no,4([[l, 2], 3], 4, 5, 6) _ 
S123S23 S123S12 
In these equations, the box coefficients &(...) are known explicitly, regardless of whether 
we are talking about MHV or NMHV amplitudes. Note that, motivated by our experience 
with the MHV six- and seven-point amplitudes, we have made the crucial assumption that 
the loop momentum dependence of the colour-dual numerators drops out after reduction to 
pentagons. Therefore, only the scalar numerators Uq remain, and the final reduction from 
pentagons to boxes is simply stated using the reduction coefficient r(. . .). We also impose 
refiection symmetry and the vanishing of all triangle numerators. 



no,6(a, b, c, d, e, /) = no,6(/, e, d, c, b, a) , 
^0,3(0, 6, c) = . 



(8.11) 
(8.12) 
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Using the Jacobi identities, we can rewrite no,5, no,4, and no,3 in terms of no,6- Eqs. (8.8)- 
(8.12) for all permutations of external legs then provide a linear system constraining the 720 
unknown quantities nQfi{a,b,c,d,e, f). This system has rank 719, subject to a number of 
consistency conditions on the box coefficients b{. . .) that are too lengthy to give here. We 
have verified that these consistency conditions are indeed satisfied by the box coefficients for 
both MHV and NMHV amplitudes, and we have solved the system for particular (generic) 
choices of kinematics. 

With 77-0,6 in hand, it is now straightforward to reconstruct all of the i dependence of 
the colour-dual numerators using the integrand oxidation procedure discussed in previous 
sections. In the MHV case, the one remaining degree of freedom is equivalent to Yuan's k 
parameter [47], and leads to a quadratic term in the hexagon numerator equal to 

71^^6(1, 2, 3, 4, 5, 6) = -3iK 6\Q) r^'' . (8.13) 

This is inconsistent with our requirement that the loop momentum dependence drops out 
after integrand reduction to pentagons, and indeed the reconstructed numerators only give 
the correct amplitude for the choice k = 0. The vanishing of n2,6 also eliminates the 
remaining degree of freedom in the NMHV amplitude. 

Seven Points 

The same setup works again at seven points. We start by writing down all of the box cut 
constraints, which relate the 5040 unknowns noj{a, b, c, d, e, /, g) to known box coefficients. 
Together with reflection symmetry and vanishing of triangle numerators, this linear system 
has rank 5019, again subject to consistency conditions on the box coefficients b{. . .). We have 
verified that the box coefficients in the MHV and NMHV amplitudes satisfy the conditions, 
and we have explicitly solved the system for 5019 of the nQ i{a, b, c, d, e, /, g) in terms of the 
21 others. 

The i dependence of the numerators can be constructed using the integrand oxidation 
procedure. As at six points, the story is not finished; the condition that i drops out of 
the numerators after reduction to pentagons provides a further 19 constraints, and we are 
left with a 2-parameter family of colour-dual numerators for seven-point MHV and NMHV 
amplitudes. For the MHV amplitude, this matches exactly the freedom found by the analysis 
in Section VII. 

We have thus established definitively the existence of colour-dual numerators for all one- 
loop six- and seven-point amplitudes in A/" = 4 SYM. Under our assumptions, the box 
numerators are unique and, for the MHV amplitude, match the box numerators obtained 
from the self-dual amplitudes. 

IX. GRAVITY AMPLITUDES 

One of the most important properties of the BCJ numerators that we constructed in 
the previous sections is the cancellation of the loop momentum dependence upon integrand 
reduction to pentagons. In this section, we will see that the same happens when we "square" 
the A/" = 4 SYM numerators to obtain the corresponding one-loop amplitudes in A/" = 8 
supergravity, following the BCJ double copy prescription, Eq. (2.7). 



27 



Six Points 



The six-point gauge theory amplitude was constructed in Section VI in the MHV sector 
and Section VIII in the NMHV sector. It is the first case where our BCJ numerators depend 
on the loop momentum. The dependence is linear, 

n6(l, 2, 3, 4, 5, 6; £) = no,6(l, 2, 3, 4, 5, 6) + ni,6(l, 2, 3, 4, 5, 6) ■ £ . (9.1) 

The dependence on £ is such that, upon integrand reduction to a pentagon, it cancels against 
that of the pentagon numerator. For instance, take the coefficient of the pentagon integral 
j^i+2)|3|4|5|6 partial amplitude; the piece proportional to ^ is 

<,(fl,2l,3,4,5,6) 

<6(1, 2, 3, 4, 5, 6) r(12|3|4|5|6) + '''^^ ' J' ' ' ' ^ ^ . (9.2) 

This cancellation can be traced back to the system of Jacobi identities that determines rii 5 
and rii 6- For a simple illustration of how tightly connected integrand reduction and Jacobi 
identities are, consider the following consistency check. Using Eq. (9.2), we get 

<6(1,2,3,4,5,6)-<6(2,1,3,4,5,6) 



<5([1,2],3,4,5,6) ^1^1 



S12 V^(12|3|4|5|6) r(21|3|4|5|6)^ 

= <5([1,2],3,4,5,6), (9.3) 

where we imposed n'l^[\l, 1], 3, 4, 5, 6) = —n'^^{\i, 2], 3, 4, 5, 6) and used the reduction coeffi- 
cients given in Eq. (5.4). 

Let us now consider the cancellation in the gravity amplitude obtained through the double 
copy, focusing again on the coefficient of the integral j^^+^^l^l^l^l^. The numerator ng is linear 
in loop momentum, so its square will have a quadratic piece, a linear piece and a scalar 
piece. The contribution of the quadratic piece is 

<6(1, 2, 3, 4, 5, 6) <,6(1, 2, 3, 4, 5, 6) r(12|3|4|5|6) 
+ <6(2, 1, 3, 4, 5, 6) <,6(2, 1, 3, 4, 5, 6) r(21|3|4|5|6) 
<5([1,2],3,4,5,6) n5'^5([l,2],3,4,5,6) 



s?2 V^(12|3|4|5|6) r(21|3|4|5|6) 

<5([1,2],3,4,5,6) nr5([l,2],3,4,5,6) 



(9.4) 

■5l2 

where we used the cancellation in gauge theory, Eq. (9.2). Notice that, because there is no 
colour ordering in gravity, we had to consider the two terms (12|3|4|5|6) and (21|3|4|5|6). 
Thus the hexagon contribution cancels the pentagon contribution. Likewise, we find a 
cancellation for the piece linear in £, 

<6(1, 2, 3, 4, 5, 6) no,6(l, 2, 3, 4, 5, 6) r(12|3|4|5|6) 
+ <6(2, 1, 3, 4, 5, 6) no,6(2, 1, 3, 4, 5, 6) r(21|3|4|5|6) 

= - ""'''^^^'^^'^'^'^'^^ K.6(l, 2, 3, 4, 5, 6) - no,6(2, 1, 3, 4, 5, 6)) 

<5([l,2],3,4,5,6)no,5([l,2],3,4,5,6) 
= . (9.5) 

S12 
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We conclude that the cancellation of ^ dependence in gravity follows from the analogous 
cancellation in gauge theory. Let us emphasise that this does not mean we can just drop 
the linear terms ni^s ■ i and ni^g ■ ^ of the numerators in all the cubic diagrams. Those terms 
ensure that the integrand reduction has the correct global properties. Consider, for instance, 
the integrand reduction in a gauge theory partial amplitude. Placing the loop momentum 
between legs 6 and 1, the hexagon integrand reduces to six distinct pentagon integrands, as 
in the expression (5.5). For five of those pentagons, the reduction proceeds as in the case 
(9.2), but for the pentagon j^^+^^l^l^l^l^^ ^^^e have 

^^ o 'i A r: a fj\ , ^5([6, 1], 2, 3, 4, 5; ^ - pe) 
^6(1,2,3,4, 5,6;£) r(61|2 3 4|5) H 

= no,6(6, 1, 2, 3, 4, 5) r(61|2|3|4|5) + ^o,5([6, 1], 2, 3, 4, 5) ^^ ^^ 

where we recall that ng is cyclically symmetric, so that 

no,6(l, 2, 3, 4, 5, 6) = no,6(6, 1, 2, 3, 4, 5) - r2i,6(l, 2, 3, 4, 5, 6) ■ (9.7) 

and 

<6(1, 2, 3, 4, 5, 6) r(61|2|3|4|5) + <5([6, 1], 2, 3, 4, 5) ^ ^ 

So we see that the naive substitution ne(l, 2, 3, 4, 5, 6; t) — > ^0,6(1, 2, 3, 4, 5, 6) does not cor- 
rectly reproduce all six pentagons simultaneously. 



Seven Points 

The situation is entirely analogous at seven points. Now the 'master numerator' riy, 
constructed in Section VII (VIII) in the MHV (NMHV) sector, is a quadratic polynomial in 
the loop momentum, 

n7(l, 2, 3, 4, 5, 6, 7; i) = ^0,7(1, 2, 3, 4, 5, 6, 7) + 711,7(1, 2, 3, 4, 5, 6, 7) ■ £ 

+ £-n2,7(l,2,3,4,5,6,7) (9.9) 

We have to consider two types of pentagons: one-mass and two-mass. 

Starting with the two-mass pentagons, and focusing on the coefficient of the integral 
j^i+2)|(3+4)|5|6|7^ cancellation of loop momentum dependence upon integrand reduction in 
gauge theory is expressed as: 

rarfi([l,2l,3,4,5,6,7) 
^^7(1,2,3,4,5,6,7) r(12|34|5|6|7) + '''^^ . J. . . . . V (l+2|34|5|6|7) 

nf6(l,2, [3,4],5,6,7) ,,,,,, , , 

+ ^'^^ ' 'I ' J' ' ' V (12|3+4|5|6|7) =0, (9.10) 



and 



S34 



<fi([l,2l,3,4,5,6,7) 
<7(1,2,3,4,5,6,7) r(12|34|5|6|7) + '''^^ . J. . . . . V (l+2|34|5|6|7) 

<g(l,2, [3,4],5,6,7) ^ , , , , ^ <5([1,2],[3,4],5,6,7) 
+ ^-^^ .U.J... r(12|3+4|5|6|7) + ^'^^^ . J. i . J. . . ; _ q ^ ^g^^^^j 

S34 S12S34 
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where we considered both the quadratic and the hnear pieces, respectively; recall that n2,5 = 
0. Since is a quadratic polynomial in £, its square is a quartic polynomial. The cancellation 
in gravity corresponds to the fact that the I dependence drops out in the following expression: 

^2(1, 2, 3, 4, 5, 6, 7; £) r(12|34|5|6|7) + ^2(2, 1, 3, 4, 5, 6, 7; £) r(21|34|5|6|7) 

+ n^(l, 2, 4, 3, 5, 6, 7; £) r(12|43|5|6|7) + ^2(2, 1, 4, 3, 5, 6, 7; £) r(21|43|5|6|7) 

n^([l,2l,3,4,5,6,7;£) , , , , , s n^rn , 21, 4, 3, 5, 6, 7; £) , , , , , ^ 
^ 6VL ' J' ' ' ' ' ' V (l+2|34|5|6|7) + '^^^ ' ^' ' ' ' ' ' V (l+2|43|5|6|7) 

■512 ■512 

^ !;I(LMMLMM,,i2|3+4|5|6|7) + !«MM^^,(2i|3+4|5|6|7) 

■534 ^534 

^ ni([l,2],[3,4],5,6,7;£) _ 

■512^534 

Notice that this requires the cancellation of pieces which are quartic, cubic, quadratic and 
linear in I. This all follows from (9.10) and (9.11), together with the Jacobi identities. 

For the other type of pentagon, namely the one-mass pentagon, let us consider the coef- 
ficient of the integral /^^+2+3)|4|5|6|7 gauge theory. The cancellations are 

nffi([l,2l,3,4,5,6,7) 
nf7(l, 2, 3, 4, 5, 6, 7) r(123|4|5|6|7) + '''^^ > J. ^ > > > r((l+2)3|4|5|6|7) 

•512 

+ ^^'^^^'^^'^^'^'^'^'^V (l(2+3)|4|5|6|7) = , (9.13) 

■534 

and 

<7(1, 2, 3, 4, 5, 6, 7) r(123|4|5|6|7) 

n''fi([l,2l,3,4,5,6,7) n't Jl, [2,31,4,5,6,7) 

+ '-'^^ ^ b . ^ . ^ V ((l+2)3|4|5|6|7) + ^''^ >[ . J> . > > V (l(2+3)|4|5|6|7) 

■512 ^534 

^ <,5([[1,2],3],4,5,6,7) ^ <5([1,[2,3]],4,5,6,7) ^ ^ 

■512^5123 ^523^5123 

For the gravity amplitude, these relations imply that the ^ dependence drops out in the 
following quantity: 

n?(l, 2, 3, 4, 5, 6, 7; £) r(123|4|5|6|7) + 71^(2, 1, 3, 4, 5, 6, 7) r(213|4|5|6|7) 

■512 .512 

+ !^l«M^^MM+eyclic(123), (9.15) 

■512^5123 

where one should add two copies of the whole expression corresponding to the other cyclic 
permutations of legs 1, 2 and 3. All the terms contribute to the coefficient of the pentagon 
integral 

X. CONCLUSION 

We have developed a systematic method of determining BCJ numerators for one-loop 
amplitudes. The strength of this method is that it makes use of the global constraints on 
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the loop momentum dependence of the numerators imposed by the Jacobi identities. In fact, 
one of the motivations for this work was the discussion by Yuan [47] of how a global constraint 
prevented the construction of six-point BCJ numerators independent of loop momentum. 
The global constraints allow us to build higher-polygon numerators starting with the box 
numerators, whence our terminology of "integrand oxidation," as opposed to reduction. 

The specific input from = 4 SYM theory that was useful is that it is possible to 
consider vanishing bubble and triangle numerators, and box numerators independent of the 
loop momentum. While we do not present a proof, our results give strong evidence that 
colour-kinematics duality holds for one-loop Af = 4 SYM amplitudes at any multiplicity and 
in any R-sector. 

We proposed a simple prescription to compute box numerators in the MHV case, using a 
very interesting connection between MHV amplitudes and self-dual gauge theory conjectured 
in [45]. We believe that this is a canonical choice of box numerators with nice properties. 
Recall that, beyond four points, it automatically prevented the 'master numerator' rim from 
depending on the highest possible power of the loop momentum, which would have been 
m — 4. The connection to the self-dual sector is suggestive, but we should also note the 
relation to the so-called MHV polygons proposed in [61]. The same set of six-point box 
numerators can be obtained by considering a simple form of MHV polygons which makes 
the cyclic symmetry manifest [47]. This particular choice and possibly other aspects of 
the Jacobi identities may also be closely related to the structure of one-loop superstring 
amplitudes studied in [8]. 

Beyond MHV amplitudes, we proposed a procedure to construct colour-dual numerators 
from knowledge of unitarity cuts alone. Using this procedure, we obtained colour-dual 
forms for six- and seven-point MHV and NMHV amplitudes in A/" = 4 SYM. In principle, 
we believe this will work for any multiplicity and in any R sector. More work is needed 
to fully understand the structure of the equations that must be solved; in particular, there 
is a set of constraints that must be satisfied by the unitarity cuts of an amplitude for a 
colour-dual form to exist. It would also be interesting to determine simple analytic forms 
for the numerators. 

The general aspects of the structure of BCJ numerators obtained through our method 
were outlined in the previous section. A very clear property of the loop momentum de- 
pendence of the numerators is that it cancels out upon integrand reduction to pentagons. 
This cancellation shows that the Jacobi identities are tightly intertwined with the integrand 
reduction coefficients. These in turn lead to cancellations in one-loop A/" = 8 supergravity 
amplitudes. While the ultraviolet properties of one-loop Af = 8 amplitudes have long been 
understood by other methods, here we have seen explicitly their intricate relationship with 
the colour-kinematics duality. It is of obvious interest to investigate the extension of this 
relationship to higher-loop amplitudes. 

As a side result of our investigations, we obtained simple trace-based formulas for integral 
reduction to pentagons and to boxes, which are more convenient in applications using the 
spinor-helicity formalism than the equivalent determinant-based formulas [57, 58]. We hope 
that these results may find a more general use. 

To conclude, our goal was to make explicit the linear map that puts a one-loop amplitude 
into a colour-kinematics dual form. We expect that the main idea explored here of using the 
global constraints on loop momentum dependence will find use beyond the one-loop level. 
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